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<N ! Abstract 

7— I . 

Recently BES2 collaboration observed an enhancement near the pp invariant 
mass spectrum. Using the covariant tensor formalism, here we provide theoretical 
^sD . formulae for the partial wave analysis (PWA) of the ip radiative decay channels 

ip — > 7pp, 7AA, 7XX, 7HH . By performing the Monte Carlo simulation, we give the 
■ angular distributions for the photon and proton in the process of J/ip — > "fpp, which 

may serve as a useful reference for the future PWA on these channels. 

o ; 

^! 1 Introduction 

J/ip and ip' radiative decay to BB (baryon and antibaryon pair) is a good channel to 
• study the possible bound or resonant states of the BB system. Abundant J/ ip and ip' 
>- ! events have been collected at the Beijing Electron Positron Collider (BEPC). More data 
will be collected at upgraded BEPC and CLEO-C. Based on the analysis of the 58 million 
J/ip events accumulated by the BES2 detector at the BEPC, recently BES2 reported 
[5] that they observed a strong, narrow enhancement near the threshold in the invariant 
mass spectrum of pp (proton - antiproton) pairs from J '/ip — > jpp radiative decays. The 
structure has the properties consistent with either an S(0 h )- or P(0 ++ )-wave Breit- 
Wigner resonance function. In the S-wave case, the peak mass is blow 2M P around 
M = 1859MeU with a total width T < 30MeV .In order to get more useful information 
about properties of the resonances such as their J PC quantum numbers, mass, width, 
production and decay rates, etc., partial wave analyses (PWA) are necessary. PWA is 
an effective method for analysing the experimental data of hadron spectrum. There are 
two types of PWA: one is based on the covariant tensor (also named Rarita-Schwinger) 
formalism [6] and the other is based on the helicity formalism [7]. Ref. [8] showed the 
connection between the covariant tensor formalism and the covariant helicity one. Ref. [9] 
provided PWA formulae in a covariant tensor formalism for ip decays to mesons, which 
have been used for a number of channels already published by BES[10, 11, 12, 13, 14, 15] 
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and are going to be used for more channels. A similar approach has been used in analyzing 
other reactions[16, 17, 18]. Ref.[19] provided explicit formulae for the angular distribution 
of photon of the ip radiative decays in the covariant tensor formalism, and also discussed 
helicity formalism of the angular distribution of the ip radiative decays to two pseudoscalar 
mesons, and its relation to the covariant tensor formalism. Now we extend the covariant 
tensor formalism to the ip — > ^BB with B represents baryons p, A, E, S etc, and derive 
the decay amplitudes for various intermediate resonant states in the framework of the 
relativistic covariant tensor formalism. 

In this paper we study the phenomenological spin parity determination of resonances. 
The plan of this article is as follows: in section 2, we present the necessary tools for the 
calculation of the tensor amplitudes, within a covariant tensor formalism. This will allow 
us to derive covariant amplitudes for all possible processes. In section 3, we present covari- 
ant tensor formalism for ip radiative decays to baryon antibaryon pairs. Since covariance 
is a general requirement of any decay amplitude, all possible amplitudes are written in 
terms of covariant tensor form. All amplitudes include a complex coupling constant and 
Blatt-Weisskopf centrifugal barriers where necessary. In section 3, we provide the angular 
distribution of the photon and proton. The conclusions are given in section 5. 

2 Prescriptions for the construction of covariant ten- 
sor amplitudes 

In this section we present the necessary tools for the construction of covariant tensor 
amplitudes. The partial wave amplitudes U- lva in the covariant Rarita-Schwinger tensor 
formalism can be constructed by using pure orbital angular momentum covariant tensors 
^w-Mi anc ^ covar i an t spin wave functions m ... Ms together with metric tensor g^ u , totally 
antisymmetric Levi-Civita tensor e^Ao- and momenta of parent particles. For a process 
a — > be, if there exists a relative orbital angualar momentum L^ c between the particle a 
and b, then the pure orbital angular momentum Lb c state can be represented by covariant 
tensor wave functions t^?)^ [7] which is built out of relative momentum. Thus here we 
give only covariant tensors that correspond to the pure S-, P-, D-, and F-wave orbital 
angular momenta: 
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where r = pb — p c is the relative four momentum of the two decay products in the parent 
particle rest frame; (f • f) = — r 2 . and 

g^{Pa) = g^u — ; (5) 

Pa 

B Lbr (Q abc ) is a Blatt-Weisskopf barrier factor [7, 21], here Q abc is the magnitude of p b or 
p c in the rest system of a, 

qLc = (Sa+ 4 s ;; Sc)2 - s b (6) 

with s a = E 2 a - p 2 . 

The spin-1 and spin-2 particles wave functions 0^(p a ,m) and (^^(pajm) satisfy the 
following conditions 

Pa<l>n{Pa, m) = 0, M (Pa, m)0* M (p a , m) = -6 mm ' 

^(P<*> m )<Pt(Pa, m) = -g^ + Pa ^ au = -g^ipa), (7) 

m Pa 

and 

P !lJ»>v>(Pa) = <l>lw(Pa, m )(t>l'u'{Pa, m) = -(g^'^u' + gav'^v^') ~ ( 8 ) 
ra Z 6 

Note that for a given decay process a — > be, the total angular momentum should be 
conserved, which means 

J a = Sfe c + L bc , (9) 

where 

S 6c = S 6 + S c . (10) 
In addition parity should also be conserved, which means 



Va = VbVc(-l) Lbc , (11) 

where r) a , i] b , and r\ c are the intrinsic parities of particles a, b, and c, respectively. From 
this relation, one knows whether L bc should be even or odd. Then from Eq. (9) one 
can find out how many different L bc — S bc combinations, which determine the number 
of independent couplings. Also note that in the construction of the covariant tensor 
amplitude, for Sb c — L bc — J a coupling, if S bc + L bc + J a is an odd number, then e^ u \ a p^ 
with p a the momentum of the parent particle is needed; otherwise it is not needed. 
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3 Covariant tensor formalism for ip decay into jBB 

The general form of the decay ip — > 7X — > jpp amplitude can be written as follows by 
using the polarization four-vectors of the initial and final states, 

A« = M^j)e:(m^ as {p b , S b ; Pc , S C )A^ = ^{mjK(™i)^s{P*, S b ; Pc , S c ) £ A^C/f^ 

(12) 

where ip^ijnj) is the polarization four vector of the ip with spin projection of mj] e v {m^) 
is the polarization four vector of the photon with spin projections of m 7 ; U- luas is the i-th 
partial wave amplitude with coupling strength determined by a complex parameter Aj. 
The spin-1 polarization vector ip^{mj) for ip with four momentum p^ satisfies 

3 

^n(. m jWv(. m j) = ~9\xv + = -g„u(p), (13) 

mj=l P 

with = 0. Here the Minkowsky metric tensor has the form 

g^ = diag(l, -1, -1, -1). 

For ip production from e + e~ annihilation, the electrons are highly relativistic, with the 
result that J z = ±1 for the ip spin projection taking the beam direction as the z-axis. 
This limits m j to 1 and 2, i.e. components along x and y. Then one has the following 
relation 

2 

J2 Vv( m </)V^'( m </) = <W(<Vi + V)- ( 14 ) 

mj=l 

For the photon polarization four vector, there is the usual Lorentz orthogonality condi- 
tions. Namely, the polarization four vector e u (m y ) of the photon with momenta q satisfies 

q v e v (m^) = 0, (15) 

which states that spin-1 wave function is orthogonal to its own momentum. The above 
relation is the same as for a massive vector meson. However, for the photon, there is an 
additional gauge invariance condition. Here we assume the Coulomb gauge in the tp rest 
system, i.e., p v e v = 0. Then we have [20] 

2 e;(m 7 ) ei ,(m 7 ) = -g^ + ^ K ^ + ^^ _ JLJL q^ = _ g (^) ( 16 ) 

with K = p — q and K u e u = 0. For X — > pp, the total spin of pp system can be either 
or 1. These two states can be represented by ip and ip a [22]. where 

ip = u( Pb ,S b )-f 5 v( Pc ,S c ), if s = 0, (17) 

T 

= u(Pb,Sb)(la )v(Pc,S c ), if S = 1. (18) 

mx + 2m 



One can see that both ip and ip a have no dependence on the direction of the momentum 
p, hence correspond to pure spin states with the total spin of and 1, respectively. Where 
Pb, p c , and Sb, S c are momenta and spin of the proton antiproton pairs, respectively. m x 
and m are the masses of X and p, p, respectively; u(pb, Sb) and v (p c , S c ) are the standard 
Dirac spinor. If we sum over the polarization, we have the two projection operators: 



J2 U a(Pb,S b )Uf3(p b ,S b ) = ( ) q 



To compute the differential cross section, we need an expression for \A\ 2 . Note that 
the square modulus of the decay amplitude, which gives the decay probability of a certain 
configuration should be independent of any particular frame, that is, a Lorentz scalar. 
Thus by using the Eqs. (14) and (16), the differential cross section for the radiative decay 
to an 3-body final state is: 

da^ 1 2 2 

"^3 z s b ,S c mj=l mi =l 



= -\ E E^j^^C.^. 

Z S b ,S c n=l 



-^EA^E^^c/r'^ E 



= E/o (20) 



where 



p 


p* 

ji 






ji 


= -\ e ^^r'^ e 



(21) 



<i<3> 3 is the standard Lorentz invariant 3-body phase space given by 



Aih ( \ x*( \ rf3q m 2 d 3 p b d 3 p c 

^3(p;g,P>,Pc) = <y (P-g-P>-Pc) (27rm(27r)3g6(27r)3ge - (22) 

4 0) = ^ (0) = 4 E vrg^xrf E « 

M=l 5 6 , 5 C 
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4m 2 ^ 

The spin sums can be performed using the completeness relations from Eq. (19): 



(23) 



Fx 



(i) 



jxva (J— L) t j*\iv' a! 



i>b + in i> c — rn 



pb + m — m 

2m ^ a 2m 



-Tr( 
i v 



mx + 2m v 2m 2m 1 mx + 2m v 2m 

r a r a i T ( $ h + m j> c -m 
(m x + 2m) 2 ^ 2m 2m 



la' ) - 

)] 



rv 



7a' ) 



b + m rf c — m 

la 



2m 



E urg^v 



(J — L) j j-*/j,u'a' 
3 



Am 2 

where we have used traces of 7 matrics, 

tr(l) = 4, tr( 7 5 ) = 0, 

tr(any odd (j 7's) = 0, 



PbaPba' + PcaPca' + PbaPca' + Pba'Pca ~ m x g aa , 



(24) 



tr{Yl u l 5 ) = 0. 



(25) 



Amplitudes for the radiative decay ip — > 7pp 

We consider the decay of a ^ state into two steps: ip — > 7 X with X — > pp. The possible 
jpc £ or x are ++ ,0~ + , 1 ++ , 2 ++ , 2 h , etc. For ^ — > 7 X, we choose two independent 
momenta p for ip and g for the photon to be contracted with spin wave functions. We 
denote the four momentum of X by K. The tensor describing the first and second steps 
will be denoted by Tf£). L and tfy...^ , respectively. 

For ip — ■> 7 ++ — > 7pp, there is one independent covariant tensor amplitude: 

U^ a = g^t {1)a . (26) 

For ip — > 7 CT + — > 7pp, there is one independent covariant tensor amplitude: 

= e^p x q a B l {Q^ x ). (27) 
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For ip — > 7l ++ — > 7J9J9, there are two independent covariant tensor amplitudes: 

Ur a = ^ vXa V\^ p „Kp§\ (28) 
ur = e ^ pxq%e a^ Kp fi) B2{Q ^ x) (2g) 

For if) — > 7l~ + , the exotic 1~ + meson cannot decay into pp. 

For ■?/> — > 72 ++ — > pp, there are three independent covariant tensor amplitudes: 

= P (2)fMUa/3 {K)tP, (30) 
Ur = g^P {2)a ^ Pm t^B 2 (Q^ x ), (31) 
^ = p(^Yp/x ] B 2 (Q^x). (32) 

For 2 ++ decaying to pp, the orbital angular momentum between the proton and antiproton 
/ could be 1 and 3; but we ignore / = 3 contribution because of the strong centrifugal 
barrier. 

For tp — > 72~ + — > pp, the possible partial wave amplitudes are the following: 

Ur = ^Pxq^B^x), (33) 
U^" = e^p xqa p,p s t^ 5 B 3 (Q^ x ), (34) 
Ut = ^ x y x q^p 5 t$B 3 (Q^ x ). (35) 

It is worth to mention here that the above partial wave amplitudes for the process 
J/t/j — > 7pp are applicable to the processes J /if) — > 7AA,7££, and 7SS as well. 



4 Monte Carlo simulation for J/ip — > jpp 

We perform Monte Carlo simulation for the decay process J/ip — ^ 7X — > 7pp, with the 
J PC of X to be ++ , (T + , 1 ++ , 2 ++ , 2 -+ . The predicted angular distributions for various 
J PC intermediate resonances could serve as a useful reference for people performing partial 
wave analysis of ip — > ^BB channels. 

BES2 reported observation of a narrow enhancement near 2m p in the invariant mass 
spectrum of pp pairs from radiative J/ip — > 7pp decays. The peak has properties consistent 
with either a J PC = CT + or ++ . We simulate these two processes with the Breit-Wigner 
mass and width of X in Ref.[5]. The angular distributions of the photon and proton are 
shown in Fig.l, while the invariant mass of pp and momentum distribution of the proton 
are shown in Fig.2. The two fits of 0~ + and ++ are indeed hardly distinguishable from 
these distributions. However, if the narrow structure is really due to a narrow CT + or ++ 
resonance, it can be distinguished by its other decay modes. While a 0~ + resonance can 
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Figure 1: Angular distributions of the photon and proton in the process J/ip — > 70 ±+ 
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Figure 2: Momentum distributions of proton and invariant mass of pp in the process 
J/ij} — > 70 ±+ — > ^ pp. The dashed curve is the case of ++ ;the solid curve is the case of 

o-+. 
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decay into rjinr, KKir, a ++ resonance cannot. On the other hand, while a ++ resonance 
can decay into two pseudoscalar mesons, such as nn, KK, a 0~ + resonance cannot. Both 
CT + and ++ resonances can decay into 4ir and rnrKK channels. Previous data on these 
channels [10, 11, 12, 13, 14, 15] have not seen such narrow structure around 2m p . This 
gives support to the explanation of pp 0~ + final state interaction [23, 24] for the observed 
narrow peak structure in pp channel only. 

The processes with other J PC intermediate X resonances are simulated by assuming 
the mass of X at 2.15 GeV with a width of 0.15 GeV. The angular distribution of photon 
and proton are shown in Figs. 3-5. The # 7 and 6 P are given in J/ip rest frame and X rest 
frame, respectively. The 9 1 angular distributions coincide with general analytic formulae 
given in Ref. [19] as it should be. One can see that various partial wave amplitudes give 
different angular distributions. By fitting the theoretical differential cross section given 
by Eq.(20) with parameters Aj to the data, one can get the magnitudes of each partial 
wave contribution. 

5 Conclusion 

In this paper we provided a theoretical formalism and a Monte Carlo study of the partial 
wave analysis for the radiative decay J/ip — * 'ypp, which are also applicable to the processes 
J ftp — > 7AA, 7SS and 7SS. We have constructed all possible covariant tensor amplitudes 
for intermediate resonant states of J < 2. For intermediate resonant states of J > 3, the 
production vertices need L > 2 and are expected to be suppressed [9]. The formulae 
here can be directly used to perform partial wave analysis of forthcoming high statistics 
data from CLEO-c and BES-III on these channels to extract useful information on the 
baryon-antibaryon interactions. 
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Figure 3: Angular distributions of the photon and proton for independent amplitudes 
of the process J/ip — > 7l ++ — > 7pp. The solid curve and dashed curve correspond to 
amplitudes (28) and (29), respectively. 
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Figure 4: Angular distributions of the photon and proton given by independent amplitudes 
Eq.(30) (stars), Eq.(31) (solid lines), and Eq.(32) (dashed lines) of the process J/ip — > 
72 ++ — > jpp. 
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Figure 5: Angular distributions of the photon and proton given by independent amplitudes 
Eq.(33) (stars), Eq.(34) (solid lines), and Eq.(35) (dashed lines) of the process J ftp — > 
72~ + — > 7pp. 
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